Genome-wide characterization of Salmonella Typhimurium genes required for the fitness under iron restriction

Background Iron is a crucial element for bacterial survival and virulence. During Salmonella infection, the host utilizes a variety of mechanisms to starve the pathogen from iron. However, Salmonella activates distinctive defense mechanisms to acquire iron and survive in iron-restricted host environments. Yet, the comprehensive set of the conditionally essential genes that underpin Salmonella survival under iron-restricted niches has not been fully explored. Results Here, we employed transposon sequencing (Tn-seq) method for high-resolution elucidation of the genes in Salmonella Typhimurium (S. Typhimurium) 14028S strain required for the growth under the in vitro conditions with four different levels of iron restriction achieved by iron chelator 2,2′-dipyridyl (Dip): mild (100 and 150 μM), moderate (250 μM) and severe iron restriction (400 μM). We found that the fitness of the mutants reduced significantly for 28 genes, suggesting the importance of these genes for the growth under iron restriction. These genes include sufABCDSE, iron transport fepD, siderophore tonB, sigma factor E ropE, phosphate transport pstAB, and zinc exporter zntA. The siderophore gene tonB was required in mild and moderate iron-restricted conditions, but it became dispensable in severe iron-restricted conditions. Remarkably, rpoE was required in moderate and severe iron restrictions, leading to complete attenuation of the mutant under these conditions. We also identified 30 genes for which the deletion of the genes resulted in increased fitness under iron-restricted conditions. Conclusions The findings broaden our knowledge of how S. Typhimurium survives in iron-deficient environments, which could be utilized for the development of new therapeutic strategies targeting the pathways vital for iron metabolism, trafficking, and scavenging. Supplementary Information The online version contains supplementary material available at 10.1186/s12863-022-01069-3.


Background
Iron is a cornerstone for numerous cellular metabolisms and serves as a cofactor for some proteins with vital functions. Iron is involved in many critical biochemical reactions, including respiration, tricarboxylic acid cycle, synthesis of metabolites, and enzyme catalysis. Therefore, iron is a crucial metal for the survival of bacterial pathogens [1,2]. The non-typhoidal intracellular S. Typhimurium can infect a wide range of hosts and cause gastroenteritis [3]. It has been estimated that S. Typhimurium is accountable for 93.8 million cases of gastroenteritis, leading to 155,000 deaths worldwide yearly [4]. As iron accessibility is vital for S. Typhimurium pathogenesis, the host uses a variety of mechanisms to sequester it from bacteria [5]. Also, it has been shown that a probiotic Escherichia coli Nissle 1917 reduces S. Typhimurium colonization by competing for iron [6]. After consuming foods or water Open Access BMC Genomic Data *Correspondence: ykwon@uark.edu contaminated with Salmonella, the pathogen reaches the intestine and breaches epithelial tissue to enter macrophages [7]. A defense mechanism that the host uses to fight against pathogens is depleting free iron via iron-sequestering proteins such as heme, hepcidin, ferritin, transferrin, and lactoferrin [8]. Hepcidin is produced in response to infection to bind and inactivate the cellular iron exporter ferroportin which causes iron concentration to decrease in the plasma and facilitates sequestration of iron in macrophages [2]. Despite a widespread counter-defensive strategy of hosts against the pathogens, S. Typhimurium thrives in the inflamed gut and can survive and replicate in macrophages [6,7]. Bacterial pathogens, including Salmonella, employ aggressive acquisition processes to scavenge iron from the hosts through the synthesis and excretion of highaffinity iron chelators named siderophores [9]. It has been also suggested that modulating host iron homeostasis may be a path to tackle multidrug-resistant intracellular bacteria [10]. Still, our understanding of the genes in S. Typhimurium that are required for survival in iron-restricted environments is incomplete.
It is highly important to characterize the entire genome of S. Typhimurium in a biologically relevant range of iron restriction to gain a comprehensive understanding of the genes and their proteins that play a role in coping with the stressor. The 2,2`-Dipyridyl (Dip) is the most commonly used, membrane-permeable iron chelator and selective agent to chelate Fe 2+ [11]. In a previous study, the promoters in S. Typhimurium that respond to 200 μM Dip were identified using a high-throughput approach based on the random promoter fusions [12]. Microarray also has been used intensively in different bacteria to profile global transcriptional responses to iron limitation, using varying concentrations of Dip: for instance, 200 μM (E. coli) [13], 160 μM (Shewanella oneidensis) [14]; 300 μM (Actinobacillus pleuropneumoniae) [15]; 40 μM (Leptospira interrogans) [16]; 200 μM (Acinetobacter baumannii) [17] and 200 μM (S. Typhimurium) [18]. RNA-seq has also been applied for transcriptomic responses to 30 μM Dip for Rhodobacter sphaeroides [19] and 200 μM Dip for S. Typhimurium [20]. In these studies, the tested bacteria were typically exposed to one selected concentration of Dip for a short time to explore the gene expression responses. On the contrary, in our recent study, we performed the selection of the genomesaturating Tn5 mutant libraries of S. Typhimurium under the iron restriction condition of varying levels of severity generated using Dip at the concentrations ranging from 100 to 400 μM [21]. The resulting Tn-seq data sets were initially analyzed to identify the genes that are essential for growth in an iron restriction-dependent manner [21]. In the current study, we re-analyzed the same data sets to identify the genes in S. Typhimurium that are required for fitness under varying levels of iron restriction.
Utilizing highly saturated Tn5 libraries and more than a quarter-billion reads from Tn5-genomic junctions, we identified the conditionally essential genes in S. Typhimurium that are required for the growth under varying levels of iron restriction. We demonstrated that sufAB-CDSE operon is important for bacterial fitness under moderate (250 μM) and severe (400 μM), but not under mild iron restriction conditions (100 and 150 μM Dip). We also found new genes that are critical for the growth under iron-restricted conditions, including the genes encoding sigma factor E and the proteins in electron transport, glycolysis and gluconeogenesis, phosphate transport, and zinc export. Finally, we also identified the genes that when deleted increase the mutant fitness under iron restriction. The genes identified in this study can be exploited as targets for the development of novel antibiotics and expand our knowledge related to iron acquisition and trafficking in S. Typhimurium.

S. Typhimurium growth response to different concentrations of 2,2`-Dipyridyl
Initially, we investigated the growth response of the wildtype S. Typhimurium 14028S to different concentrations of iron chelator Dip. The examined Dip concentrations ranged from 100 to 2000 μM. As illustrated in Fig. S1, the final optical density (OD 600 ) of the bacterial cultures after 18 hr. incubation at 37 °C reduced as the concentration of Dip increased. The bacteria did grow in the presence of Dip at the concentrations of 100 to 500 μM. But at 1000 μM Dip and above the bacteria could hardly grow with only a marginal increase in the optical density at 1000 μM. We found a significant decrease of OD 600 in the presence of 100 μM Dip as compared to the control culture with no Dip (p < 0.05). This was an indicator that the 100 μM Dip had a negative effect on the growth of S. Typhimurium. Therefore, we decided to use Dip concentrations ranging from 100 to 400 μM for the following Tn-seq selections; the concentrations of 100, 150, 250, and 400 μM chosen for selections of the Tn5 library are hereafter referred to Dip100, Dip150, Dip250 (Dip250-I and Dip250-II for 2 independent Tn-seq selections), and Dip400, respectively (Fig. S2). As the concentration of Dip increased, the growth rate reduced, and maximum OD 600 decreased ( Fig. 1A and Fig. 1B, respectively). The Dip250 had a profound effect on the growth and it reduced maximum OD 600 by 34% in reference to LB control (Table S1). Dip400 had a severe effect on S. Typhimurium growth; it reduced the growth rate by 26% and maximum OD 600 by 48% as compared to the control. Since Salmonella encounters host niches with different concentrations of available iron, we reasoned that our Tn-seq selection conditions representing a wide range of iron restrictions are more relevant in revealing the strategies Salmonella employs to cope with iron-restriction stress during infection in the host as compared to one condition with a fixed level of iron restriction.

Summary of Tn-seq reads
The study design for the mutant selections is illustrated in Fig. S2. The respective controls were LB-II, and LB-III and the varying levels of iron-limited conditions were Dip100, Dip150, Dip250-I, Dip250-II, and Dip400. LB-II and LB-III required about 5 hr. to reach the mid-exponential phase, while this time was extended as the concentration of Dip increased. It took 10 hr. and 12 hr. for Dip250 and Dip400, respectively (Table S2). We successfully mapped more than 173 million sequence reads to the chromosome of S. Typhimurium 14028S (NC_016856.1) for all selected conditions combined. The mean length of mapped genomic junction sequences was 91 nucleotides long. The highly saturated mutant library (Library-AB) contained 193,728 unique insertions (Table S3). The high number and long length of the mapped reads allowed us to identify conditionally essential genes with high precision. Previously, we showed that our Tn-seq protocol is highly reproducible [22]. A comparison of the Tn-seq profiles between the two biological replicates, Dip250-I and Dip250-II, indicated that the correlation coefficient (r) of unique insertions per ORF was 0.995 (Fig. S3A). In addition, the r of essentiality indices per ORF as calculated by Tn-seq Explorer [23] between Dip250-I and Dip250-II was 0.990 (Fig. S3B). These results indicate the robustness and reproducibility of our Tn-seq library protocol used in this study.

Genes implicated in the growth under iron-restricted conditions
Combining the results from all Tn-seq analyses under varying levels of iron-restricted conditions, we identified 58 genes for which the deletion mutants displayed either increased or reduced fitness in response to different concentrations of iron chelator Dip. Mutant fitness was reduced for 28 genes, while increased for 30 genes (Supplementary Data set 1). In other words, the 28 conditionally essential genes are required for the robust growth of S. Typhimurium under iron-restricted conditions. On the contrary, the 30 genes exhibited increased mutant fitness when deleted. We also identified essential genes of S. Typhimurium under these conditions, and the details of the findings were reported previously [24].

Conditionally essential genes required for fitness under iron-restricted conditions
Our results of genome-wide analyses using Tn-seq show that 28 conditionally essential genes of S. Typhimurium are required for growth under different levels of iron restriction. Some of the genes have been implicated to have a role in growth under iron-restricted conditions in previous studies, but other genes have never been associated with iron restriction (Table 1).
Iron-sulfur cluster assembly genes. We found that iron sulfur cluster assembly operon sufABCDSE is required for the growth of S. Typhimurium under ironrestricted conditions ( Table 1). The fitness of all six genes in suf operon was reduced significantly. Although the genes in suf operon are required for the bacteria to combat the iron-deficient environments, the fact that numerous reads corresponding to the Tn5 insertions in these genes were detected indicates these genes are not  Table 1 The genes in S. Typhimurium 14-028S that are required for fitness under iron restriction conditions Dip represents 2,2`-Dipyridyl and the numbers indicate the concentrations of Dip (μM). The mutants with significantly reduced fitness are shaded with green color (p < 0.05). In the RNA-seq study (reference 20), S. Typhimurium was exposed to 200 μM Dip for 10 min in Lennox broth Karash et al. BMC Genomic Data (2022) 23:55 E. coli and Salmonella have two Fe-S cluster assembly systems, isc, and suf. Under iron-restricted conditions, E. coli utilizes suf operon for Fe-S cluster assembly [21]. sufC of Salmonella Typhi was previously shown to be required for survival in macrophages [25]. The ortholog of sufS in Mycobacterium tuberculosis is implicated in iron metabolism because ΔsufS mutant of M. tuberculosis exhibited longer doubling time in the presence of Dip [26]. In addition to the vital role of suf system, our Tn-seq data also indicated that iscR is required for S. Typhimurium growth under iron-restricted conditions ( Table 1). The suf operon is controlled by iscR in E. coli [27]. IscR is not only involved in Fe-S cluster biogenesis but also implicated as a pleiotropic transcriptional regulator. S. Typhimurium iscR regulates SPI-I TTSS genes (Salmonella pathogenicity island-1 type III secretion system) when the iron level is low in vivo [28]. Pseudomonas aeruginosa IscR is considered a global regulator and it senses Fe-S cluster proteins. P. aeruginosa possesses only iron-sulfur cluster (ISC) system [29] and ΔiscR mutant showed attenuated virulence in Drosophila melanogaster and mouse peritonitis models [30]. Further, it was found that IscR regulates the expression of more than 40 genes involved in Fe-S cluster homeostasis in E. coli [31]. Therefore, we speculate that IscR is required for fitness in S. Typhimurium through its regulation of suf operon. We conducted a protein-protein interaction network analysis and the result also indicates that four Suf proteins (SufB-CDS) interact with IscR (Fig. S4). The molecular mechanisms of Fe-S proteins in bacteria have been extensively reviewed elsewhere [32]. Moreover, nfuA (yhgI) encodes an Fe-S cluster carrier protein, which was reported as a scaffold/chaperone for damaged Fe-S cluster proteins in E. coli [33]. Our Tn-seq screening shows that nfuA is required for S. Typhimurium growth under the ironrestricted condition. The ΔnfuA mutant in P. aeruginosa was sensitive to Dip and less virulent in C. elegans [34]. Deletion of nfuA in E. coli also caused the mutant to be susceptible to iron depletion [35]. The ortholog of nfuA in Acinetobacter baumannii plays a role in intracellular iron hemostasis and the bacterium cannot grow in ironchelated media when the gene is deleted [36]. Here we demonstrate for the first time that nfuA is also critical for S. Typhimurium growth under iron-restricted conditions, and we speculate the protein is possibly involved in Fe-S cluster biogenesis. Iron homeostasis. NAD(P)H-flavin reductase, fre, was also required under iron-restricted conditions (Table 1). Fre protein likely reduces free flavins, and consequently, the lower level of flavins reduces ferric iron to ferrous iron in E. coli [37,38]. We speculate that the protein product of fre gene does the same function in S. Typhimurium by reducing the Fe + 3 of siderophores to F + 2 from fepBDGC system. The ggt is another gene that S. Typhimurium uses to cope with iron restriction. The product encoded by ggt (γ-glutamyltranspeptidase) is an important enzyme in glutathione metabolism and is required for fitness under iron-restricted conditions. It has been suggested that ggt plays a role in Fe-S cluster biosynthesis in Saccharomyces cerevisiae [39]. In Campylobacter jejuni, ggt contributes to the colonization of gut in chicken and humans [40]. Helicobacter pylori lacking ggt fails to persistently colonize the stomach of mice [41]. The role of glutathione in iron trafficking was previously reported in S. Typhimurium and E. coli [42,43] but to the best of our knowledge this is the first report on the ggt role in Salmonella growth in iron restricted conditions. We also identified the genes that import iron from extracellular media. A siderophore gene fepD encodes iron enterobactin transporter membrane protein in E. coli [44]. fepD was required for S. Typhimurium growth in iron-restricted conditions. In S. Typhimurium FepD Fig. 2 sufS is required for the optimal growth of S. Typhimurium under iron-restricted conditions. S. Typhimurium 14028S wild-type and ΔsufS mutant were grown in LB broth, supplemented with 0 (control), 100 or 150 μM 2,2`-Dipyridyl (Dip). OD 600 were recorded every 15 minutes during incubation at 37 °C for 24 hr. in a 96-well plate to show the growth responses (A) and to calculate the doubling times (B). Data represent at least three replicates. Statistical significance of doubling times was determined by unpaired two-tailed t test, ***p < 0.001 is a part of FepDGC ABC transporter and is involved in the uptake of siderophore salmochelin [45]. Previously, we showed that fepD is important for the bacterium to resist oxidative stress [22]. FepD interacts with TonB and ExbB (Fig. S4). tonB and exbB were also required for S. Typhimurium growth under iron-restricted conditions ( Table 1). It has been suggested that siderophore complexes depend on TonB to energize the active transport across the membrane via TonB-ExbB-ExbD complex [46]. In S. Typhimurium tonB-mediated iron uptake is involved in the colonization of this pathogen in the Peyer's patches and mesenteric lymph nodes of mice [47]. This complex interacts with Suf system via interactions of SodA-NfuA (Fig. S4).
Sigma E factor. We found ropE and degS are required for S. Typhimurium fitness under iron-restricted conditions. rpoE encodes RNA polymerase sigma E factor, while degS encodes serine endoprotease. In E. coli, rpoE and degS are essential genes; RpoE is an extra-cytoplasmic factor that activates in response to envelope stress. The activation starts by unfolding outer membrane proteins (OMPs) and ends with proteolysis of anti-sigma E factor by DegS to free RpoE and initiate transcription [48,49]. In S. Typhimurium ropE and degS are not essential genes [24]. RpoE responds to a variety of extra-cytoplasmic stresses in bacteria and the role of this sigma factor has been determined for pathogenesis and virulence; in Salmonella the expression of rpoE is activated by different types of stressors [50,51]. Remarkably, our findings demonstrated that rpoE insertion mutant is attenuated completely under severe iron-restricted conditions (Dip400). This is reflected in the fact that no sequence reads were detected under Dip 400 (Supplementary Data set 2). To confirm the phenotype of these two genes, we grew single deletion mutants of degS and rpoE in LB supplemented with varying concentrations of Dip. After 24 h growth in a 96-well plate, we did not observe a significant change in growth rate and Max OD 600 between the mutants and wild-type. Therefore, we used spot dilution assay to confirm Tn-seq results for degS and rpoE in S. Typhimurium. As expected, S. Typhimurium degS and rpoE mutants exhibited growth defects in the presence of Dips as compared to the wild-type (Fig. 3). Whereas Dip100 did not exhibit the expected phenotype, ΔdegS and ΔrpoE formed smaller colonies at Dip200 as compared to the wild-type. At Dip300, both mutants showed a strong phenotype as compared to the wild-type. At Dip400, ΔdegS and ΔrpoE did not grow while the wildtype did grow slowly. The results of this spot dilution assay confirm that degS and rpoE are playing important roles in S. Typhimurium growth under iron-restricted conditions. The lack of discernable difference in the growth phenotype in broth media, which was captured in the spot dilution assay, indicates the high resolution of our Tn-seq assay in detecting small differences in the mutant fitness within the population of Tn5 mutants. It was previously shown that S. Typhimurium lacking degS survives very poorly in the macrophages and was slightly attenuated in mice as compared to the wild-type strain, whereas S. Typhimurium ΔrpoE mutant was attenuated more than 500-fold as compared to degS mutant in the mice [50,52]. In addition, degS plays an important role in S. Typhimurium survival in elevated temperatures and is required for full virulence [53]. rpoE can be activated by acid stress in S. Typhimurium and the gene contributes to bacterial survival in the acidified phagosomal vacuole [54]. Microarray analysis indicates that 58% of S. Typhimurium genes are affected by rpoE and there is a strong connection between SPI-2 and rpoE [55]. Also, it has been proposed that Salmonella Typhi invasion and intracellular survival are underpinned by rpoE via regulation of SPI-1 and SPI-2 [56]. Moreover, it has been suggested that rpoE regulates antibiotic resistance in S. Typhi through downregulation of the OMP genes and upregulation of the efflux system [57]. Lastly, pertussis toxin (PT) and adenylate cyclase toxin (ACT) are arsenals utilized by Bordetella pertussis to kill and modulate host cells and the expression of these two toxins is indirectly modulated by rpoE [58]. Altogether, rpoE has a broad impact on bacterial fitness and survival in the presence of various host stressors. To the best of our knowledge, this is the first report on the role of rpoE in S. Typhimurium during iron starvation.
Other miscellaneous pathways. In addition to the genes described in the previous sections, our Tn-seq analysis also identified numerous genes important for fitness under iron-restriction conditions, which are associated with other various biological pathways. These other pathways include (1) oxidative stress, (2) porphyrin biosynthesis, (3) electron transport, (4) gluconeogenesis, (5) glycolysis, (6) osmotic stress, (7) phosphate transport, and (8) Zinc exporter.
Our analysis indicates that sodA was required for S. Typhimurium growth under iron-restricted conditions. sodA gene encodes superoxide dismutase which detoxifies reactive oxygen species. In E coli the gene is under the control of Ferric Uptake Regulation (Fur) [59]. S. Typhimurium ΔsodA mutant showed a reduced capacity to invade HeLa cells and form biofilm as well as to resist chicken serum and reactive oxygen species [60]. It was also previously shown that iron restriction can induce sodA expression in vitro in S. Typhimurium [47]. hemL encoding Glutamate-1-semialdehyde 2,1-aminomutase is required for S. Typhimurium growth under ironrestricted conditions. The gene is involved in the biosynthesis of 5-aminolevulinic acid from glutamate via the five-carbon pathway [61]. ydgM (rsxB) is also required for S. Typhimurium growth under iron-restricted conditions (Table 1). It has been suggested that the product of rsxB plays a major role as the core electron mediators to reduce SoxR (redox-sensitive transcriptional activator), and it is a part of a reductase complex located in the cytoplasmic membrane in E. coli [62]. The gene ndh encoding type II NADH:quinone oxidoreductase (NAHD-II) is required for S. Typhimurium to cope with iron-restricted conditions. NAHD-II is a membrane-bound dehydrogenase that plays a central role in respiratory chains in many prokaryotes [63]. Our analysis shows that Ndh interacts with GpmA (Fig. S4). Fitness of ΔgpmA was reduced under iron-restricted conditions. The gene encodes a phosphoglycerate mutase which is involved in glycolysis and gluconeogenesis. These two genes (ndh and gpmA) belong to different pathways, and how they interact under iron-deficient environments warrants further investigations. It has been suggested that OsmE is a putative protein regulated by osmotic stress in E. coli [64]. Our Tn-seq shows that osmE is required for S. Typhimurium growth under iron-restricted conditions. pstAB is also required for S. Typhimurium growth under iron-restricted conditions. These genes are parts of ABC transporter pstSACB complex which contributes to phosphate import. pstB provides energy to the phosphate transporter system via ATP hydrolysis [65]. In E. coli, pstSACB is upregulated in response to phosphate-limited conditions [66]. It is unclear currently how pstAB is important for fitness under iron-limited conditions. zntA is required for S. Typhimurium to grow in iron-restricted conditions. The gene encodes a zinc exporter, which has a role in zinc homeostasis [67]. zntA has been implicated in the resistance of S. Typhimurium to zinc and copper and is critical for its full virulence [68]. Recent work showed that S. Typhimurium ZntA is a part of zinc efflux transporter required to diminish cytotoxic effects of free zinc and to resist nitrosative stress [69]. It was unexpected that our result showed zntA is critical for the fitness of this pathogen in iron-deficient media. To validate this finding, S. Typhimurium lacking zntA was grown in LB media supplemented with 100 and 150 μM Dip (Dip100 and Dip150). The results confirmed our Tn-seq finding and showed that zntA is important for the growth of S. Typhimurium in vitro because the growth rate of the mutants was reduced under Dip stress in comparison to the control without Dip (Fig. 4). Our Tn-seq results indicate that zntA is required for all tested conditions, from the mild to the severe iron restriction (Table 1). Regarding the specificity of the used Dip in this study, it has been suggested that Dip cannot chelate zinc, excluding the possibility that the requirement of zntA is through depletion of zinc [70].
Uncharacterized genes. We also found four previously unknown genes to be important for the growth of S. Typhimurium under iron-restricted conditions: STM14_4330, STM14_4612, ygjQ, and yhfK. These genes are annotated in UniProt databases as follows. STM14_4330 putative sugar kinase; STM14_4612 putative cytochrome c peroxidase; ygjQ putative integral membrane protein; yhfK putative inner membrane protein.

Dynamics of the genetic requirements in response to iron restriction levels
In this work, we exposed S. Typhimurium Tn5 libraries to different concentrations of iron chelator Dip for comprehensive identification of all genes that contribute to the fitness of this pathogen under iron restriction. For convenience, we categorized the 4 iron-restricted conditions we used into 3 levels: mild (Dip100 and Dip150), moderate (Dip250), and severe iron-restriction (Dip400). Since S. Typhimurium faces iron limitation at different levels of severity during infection in the host, we reason that our strategy for Tn-seq selection under a wide range of iron restrictions can comprehensively capture the genes that are important in coping with the iron stressor at various stages or niches of Salmonella infection in the host. Interestingly, we found that not all 28 conditionally essential genes required for fitness under iron restriction were required under all 4 levels of ironrestricted conditions ( Table 1). The 28 genes were clustered into different groups according to the patterns of the mutant fitness values in response to different levels of iron restriction (Fig. S5 and Fig. S6). The genes of Suf system, sufABCDSE, which were required under moderate and severe conditions (Dip250 and Dip400), were not required for mild iron restriction conditions (Dip100 and Dip150). The fitness of these mutants was reduced more when iron-restriction severity elevated: the average fitness of these mutants was − 1.8 and − 2.6 for Dip250 and Dip400, respectively. This suggests that S. Typhimurium uses Suf system to survive in moderate and severe iron restricted niches., The fitness of the mutant in iscR encoding sufABCDSE regulator was reduced in Dip250 and Dip400, − 2.75. However, S. Typhimurium iscR insertion mutants grew better in mild iron-restricted conditions (Dip100 and Dip150) with a fitness of 3.09. For the nfuA, the gene was required only under moderate ironrestricted conditions.
Interestingly, the fitness of the siderophore gene fepD was − 2.91 in Dip100 and reduced to − 4.27 in Dip150 and − 3.56 in Dip250, while fepD was not required in D400 for S. Typhimurium growth. This suggests that fepD is critical for bacterial growth under mild and moderate iron-restricted conditions but not under severe iron restriction. The tonB insertion mutant exhibited similar phenotype as fepD with fitness of − 6.43, − 4.32 and + 1.44 for Dip150, Dip250 and Dip400, respectively. exbB insertion mutant also behaved similarly. Unexpectedly, the sigma E factor had the highest reduced fitness in moderate iron-restricted conditions (the fitness − 5.97) and it reduced further in severe iron-restricted conditions (− 7.24). Also, the fitness of degS insertion mutant, the mutant without functional anti-sigma E factor, was − 2.60 for Dip250, and degS insertion mutant recorded the lowest fitness of − 8.35 in Dip400. This demonstrates that rpoE and degS are increasingly required for S. Typhimurium growth under both iron-restricted conditions in a manner dependent on the concentration of Dip. The 6 genes, sodA, gpmA, ydgM, yhfK, STM14_0026, and STM14_4612, were only required in moderate iron-restricted conditions for S. Typhimurium growth. Fig. 4 zntA is required for the growth of S. Typhimurium under iron-restricted conditions. S. Typhimurium 14028S wild-type and ΔzntA mutant were grown in LB broth supplemented with 0 (control), 100 or 150 μM 2,2`-Dipyridyl (Dip). OD 600 was recorded every 15 minutes during incubation at 37 °C for 24 h in a 96-well plate. Statistical significance was determined by unpaired two-tailed t test, *P < 0.05 Whereas ggt, hemL, ndh, osmE, STM14_4330, and ygjQ were only required in severe iron-restricted conditions. The fitness of pstB in moderate iron-restricted conditions was − 1.34 and in severe iron restriction, it became − 2.04 while pstA was required only in severe iron restriction conditions for S. Typhimurium growth. Finally, zntA was significantly required in Dip150 and Dip400, with the fitness values of − 1.39 and − 1.03. Altogether, these findings demonstrate the requirements of the genes identified in this study for fitness are dependent on the severity of iron restriction and suggest that certain genes might be required for fitness under a specific range of iron restriction.

Genes that increase mutant fitness under iron-restricted conditions upon inactivation
In this study, our main goal was to identify the genes that when deleted reduce the fitness under iron-restricted conditions. However, we also identified 30 genes that when deleted result in increased fitness, and their mutants grow better under iron-restricted conditions. These phenotypes were mainly observed in moderate (Dip250) and severe iron-restricted conditions (Dip400) (Supplementary Data set 1 and Table 2). We briefly categorize them based on their biological functions and highlight a few of these genes. First, the genes involved in nucleotide biosynthesis and metabolism: the fitness of the subunits of DNA polymerase V, umuD was increased under iron-restricted conditions. DNA polymerase V incorporates 8-oxo-guanine into DNA during replication, it has been shown that E. coli lacking umuD confers resistance to the antibiotics and can grow in the presence of the antibiotics [71]. In addition, the fitness of guaB and purA insertion mutants also increased. These two genes catalyze the first step in the de novo synthesis of guanine and adenine from inosine 5′-phosphate (IMP). Second, the genes involved in TCA cycle: we found that the mutant fitness of acnB, icdA, and sdhD genes, which encode TCA cycle enzymes, increased under iron-restricted conditions; S. Typhimurium lacking one of these genes can grow better under iron-restricted conditions, whereas the fitness of other mutants in TCA cycle did not change. A similar occurrence was observed when E. coli was exposed to lethal doses of antibiotics. The survival of the deletion mutants in acnB or icdA increased against bactericidal antibiotics [11]. Third, the genes involved in carbohydrate metabolism: S. Typhimurium strains with a deletion in yaeD, eda, or STM14_2709 can grow better under iron-restricted conditions. Eda encodes Entner-Doudoroff aldolase which has a central role in sugar acid metabolism and detoxification of metabolites in E. coli [72]. Finally, the genes involved in various pathways or functions, including lipopolysaccharide biosynthesis (rfbB, rfaC, rfbH, and rfbI), integral component of membrane (eyeiB, ychH, STM14_0726), membrane transporters (sapA, ompW, smvA), membrane-localized protease (htpX) and outer membrane protein assembly factor (nlpB). This suggests that mutations in any of the twelve genes associated with the membrane-related functions led S. Typhimurium to grow better under iron-restricted conditions. Among all those genes, only three genes (acnB, rfbI, and STM14_0726) possibly bind to Fe-S. Currently, we are lacking a scientific explanation of how the deletions in these genes caused the bacteria to grow better under iron-restriction stress.

Conclusions
In this work, we characterized the genome of S. Typhimurium at a system-wide level to identify the genes required for growth under the stress of iron restrictions. We also demonstrated the requirements of these genes for fitness alter according to the severity of the iron restriction. We validated the phenotypes only for 4 single deletion mutants for sufS, zntA, degs, and rpoE genes. However, our previous works indicated that 84% (50 mutants tested) of the genes identified by our Tn-seq method and bioinformatic pipeline displayed the expected phenotypes. We plan to move with these in vitro findings for further evaluation using macrophage cell lines and mice. The results of this study can be exploited for the development of effective therapeutic strategies and it can expand our knowledge about how Salmonella survives in ironrestricted environments.

Growth response of S. Typhimurium to different concentrations of 2,2`-Dipyridyl (dip)
A single colony of Salmonella enterica serovar Typhimurium ATCC 14028S was inoculated into a two ml LB (Luria-Bertani) broth medium in a five ml tube and incubated overnight in dark for 16 h. The single colony was taken always from an agar plate with the colonies not older than 10 days. The next day, freshly prepared LB broth media supplemented with different concentrations of Dip were inoculated with the S. Typhimurium culture at 1:200 dilutions. Dip was dissolved in ethanol and then diluted in autoclaved MQ-H 2 O before adding it to LB broth media. The cultures (200 μl/well) were immediately added to a 96-well microplate, and incubated in a Tecan Infinite 200 microplate reader at 37 °C with a shaking amplitude of 1.5 mm and shaking duration of five s. OD 600 of the cultures were measured every 10 min. After 18 h incubation, the data were collected and run on GrowthRates script to calculate the growth rate, and maximum OD 600 [73].

Preparation of Salmonella Typhimurium Tn5 mutant libraries
The Tn5 mutant libraries were constructed as previously described [22,74]. Briefly, a spontaneous mutant strain of S. Typhimurium ATCC 14028S, which is resistant to nalidixic acid (NA R ), was mutagenized by conjugation utilizing Escherichia coli SM10 λpir harboring a suicide pBAM1 transposon-delivery plasmid vector (Amp R ) as Table 2 The genes in S. Typhimurium 14,028 s that upon deletion increased the fitness of the mutants under iron restriction conditions Dip represents 2,2`-Dipyridyl and the numbers indicate the concentrations of Dip (μM). The mutants with significantly increased fitness are shaded with green color (p < 0.05) the donor cell [75]. Overnight growth cultures from the donor and recipient bacteria each were mixed, concentrated on the nitrocellulose filter, and incubated for five h at 37 °C. The conjugants were plated, then colonies scrapped, and stored.

Tn-seq selection conditions of iron restrictions
We constructed a transposon Tn5 mutant library which consists of 325,000 random mutants as described previously [22]. These transposon mutants were recovered on 50 LB agar plates (Library-A). After sequencing, it turned out that there was an insertion in every 42 nucleotides on average in the S. Typhimurium chromosome and 90% of ORFs had insertions. Then, we made another library containing 325,000 Tn5 mutants (Library-B) and combined it with Library-A, forming Library-AB. As a result, there was an insertion per 25 nucleotides on average and 92.6% of ORFs had insertions in Library-AB.
The aliquot of the Tn5 mutant library was thawed at room temperature on ice and then diluted in LB broth. The library was incubated at 37 °C with shaking at 225 rpm for an hour and then washed twice with PBS. This allows the mutants to reactivate and get rid of DMSO residues. The library was inoculated to 20 ml LB broth in a 300 ml flask and LB broth supplemented with 100, 150, 250, or 400 μM Dip (Dip100, Dip150, Dip250, and Dip400, respectively), and LB broth without Dip was used as a control. The inoculum of each library represented approximately 10 cells for each mutant. Library-A was used for selection conditions Dip100, Dip150, and LB-II (Control), while Library-AB was used for the conditions of Dip250-I, Dip250-II, Dip400, and LB-III (Control). The cultures were incubated at 37 °C with shaking at 225 rpm in a dark and humidity-controlled incubator till they reached the mid-log phase. Then, the cultures were pelleted and stored at − 20 °C. The genomic DNA was extracted from each culture, and Tn-seq amplicon libraries were prepared for HiSeq sequencing as described previously [74].

Tn-seq library generation for HiSeq sequencing
Previously we developed a robust method for Tn-seq library preparation, and therefore the Tn-seq amplicon libraries for HiSeq sequencing were prepared according to the protocol [22,74]. Briefly, genomic DNA was extracted and subjected to a linear PCR to fish out genomic junctions, utilizing Tn5-DPO primer (Table S4). Then, amplified linear junctions C-tailed. Next, the exponential PCR was conducted utilizing P5-BRX-TN5-MEO and P7-16G primers (Table S4). Finally, size selection of the library was conducted on agarose gels, libraries were sequenced on a HiSeq Illumina platform.

Tn-seq data analysis
The output data of the Hi-Seq sequencer were downloaded onto the High-Performance Computing Center (AHPCC) at the University of Arkansas. Since samples were multiplexed before sequencing, a custom Python script was used for de-multiplexing them. The script searched for the six-nucleotide barcode of each library with no mismatch allowed. Tn-Seq Pre-Processor (TPP) tool was utilized to extract transposon genomic junctions [76]. We modified the script of the TPP to process our sequences. In a fixed sequence window, TPP searched for the 19-nucleotide Tn5 inverted repeat (IR) and identified the five nucleotides GACAG at the end of the IR sequence. The genomic junctions of S. Typhimurium that start immediately after GACAG were extracted and the C-tails were trimmed. The genomic junction sequences of less than 20 nucleotides were excluded and the remaining junction sequences were mapped to Salmonella enterica serovar Typhimurium 14028S genome utilizing BWA-0.7.12 [77]. Essentiality Indices (EI) were calculated using Tn-seq Explorer [23].

Identification of conditionally essential genes required for growth under iron-restricted conditions
The conditionally essential genes were identified utilizing TRANSIT tool [76]. The resampling algorithm was utilized for the analyses in TRANSIT. The LB-II and LB-III were the inputs (controls), while Dip100, Dip150, Dip250-I, Dip250-II, and Dip400 were the outputs (experiments). Trimmed Total Reads (TTR) were used as a normalization method and 10,000 samples were used for the analysis. The insertions in 5% of N-terminal and 10% of C-terminal of ORFs were excluded. Genes were considered conditionally essential when p values were < 0.05.
Phenotypic analyses of the single deletion mutants. For assessment of the growth phenotypes of the mutants in LB broth, the overnight cultures of the wild-type S. Typhimurium, ∆sufS, and ∆zntA were inoculated into LB broth containing 0, 100, or 150 μM Dip with the inoculum diluted at 1:200. Then, 200 μl of cultures were directly added into 96-well microplates and incubated in Tecan infinite 200. The incubation time was 18 h at 37 °C. The OD 600 data was collected every 10 min and used to calculate the growth rate with GrowthRates script [73] and to obtain the maximum OD 600 . For the spot dilution assay, the overnight cultures of the wild-type, ∆degS, and ∆rpoE were serially diluted from 10 0 to 10 − 7 in a 96-well plate. LB agar plates were prepared 2 days earlier to contain 0, 100, 200, 300, and 400 μM Dip. Five μl of serially diluted cultures were spotted on the agar plates and let dry completely at room temperature. The plates were incubated at 37 °C and results were recorded after 24 h. All single deletion mutants were obtained through